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Most of the current theories on the p-wave superfluid in cold atomic gases are based on the 
effective-range theory for the two-body scattering, where the low energy p-wave scattering amplitude 
fi{k) is given by fi(k) = —l/[ik + l/(Vk 2 ) + 1/1Z], where k is the incident momentum, and V 
and 7Z are the fc-independent scattering volume and effective-range, respectively. However, due 
to the long-range nature of the van der Waals interaction between two colliding ultracold atoms, 
thep-wave scattering amplitude of the two atoms is not described by the effective-range theory 
[J, 0- In this paper we provide an explicit calculation for the p-wave scattering of two ultracold 
atoms near the p-wave magnetic Feshbach resonance (PMFR). We show that the low energy p- 
wave scattering amplitude in the presence of PMFR takes the form fi(k) = —l/[ik + l/(V eff fe 2 ) + 
l/(S eS k) + l/1Z efi ] where V ctt , <S cff and lZ eS are fc-dependent parameters. Based on this result, we 
show sufficient conditions for the effective range theory to be a good approximation of the exact 
scattering amplitude. Using these conditions we show that the effective-range theory is a good 
approximation for the p-wave scattering in the ultracold gases of 6 Li and 40 K when the scattering 
volume is enhanced by the resonance. 

PACS numbers: 



I. INTRODUCTION 

Recently, ultracold atomic gases with strong p-wave in- 
teraction have attracted broad interest both experimen- 
tally [Ml and theoretically [HH. The 

p-wave mag- 
netic Feshbach resonances, which can generate tunable 
p-wave ineratomic interactions, have been observed in 
the cold gases of 40 K Sjl, 6 Li 6 Li- 87 Rb mixture 

[H[, 6 Li- 40 K mixture 13, and 40 K- 87 Rb mixture 
The p-wave Feshbach molecules have also been created 
and studied in the gases of 40 K @, 0] and 6 Li @, M~ 
OjJ • These experimental achievements stimulate theoret- 
ical researches on the quantum superfluid in ultracold 
atomic gases with strong p-wave interactio ns |16h 52||. as 



well as the relevant few-body problems [52--6E 

Until now, most theories of ultracold atomic gases 
with strong p-wave interactions (35l - [38l [iol - li^ . [52rl57j are 
based on the low energy expansion of the p-wave scatter- 
ing amplitude /i(fe) given by the effective-range theory 



fi(k) 



ik + vk 2 



i_ 

■R 



(1) 



Here k is the relative momentum of the two atoms; V 
is the scattering volume and TZ is the effective-range. 
The effective-range theory for the scattering amplitude 
is used in both theories of p-wave atomic superfluids 
J3B-S, I4(il - l43 . [52| and related few-body problems l52h - 
l57j . In particular, the separable two-body potential [37| 
V(k, k') — Xiw(k)w(k / ) used in the many-body Hamilto- 
nian for the p-wave atomic superfluid is derived directly 
from low energy expansion {1} of the scattering ampli- 
tude. 



However, the effective-range theory is correct only for 
the short range potentials (e.g., Yukawa potential) 
which decays faster than any power function r -7 in the 
large interatomic distance limit r — > 00. Here r is the 
relative coordinate between two atoms. For a realistic in- 
teraction between two cold atoms, which is described by 
a long-range potential dominated by the van der Waals 
term — fi, 2 /3|/(r 6 m) in the limit r — > 00. Here m is the 
single-atom mass, H is the Plank constant and Pq is the 
van der Waals length. Due to the long-range van der 
Waals potential, the effective-range theory and the low- 
energy expansion (fTJ) of the p-wave scattering amplitude 
is not applicable any longer [l|, 0] . 

In the presence of a p-wave magnetic Feshbach res- 
onance (PMFR) in the ultracold gases of polarized 
fermionic atoms, the p-wave scattering amplitude of the 
atomic collision is contributed by both the background 
potential in the open channel and the bound state in the 
closed channel. The long-range nature of the background 
potential makes the final scattering amplitude to be in- 
consistent with the effective-range theory. 

Therefore, it is essential to investigate the condition 
under which the effective-range theory (TTJ) can be used 
as an approximation of the exact p-wave scattering am- 
plitude under a PMFR. If the effective-range theory pro- 
vides a good approximation of the scattering amplitude, 
then the previous theories on p-wave superfluid would 
be applicable; nevertheless, if the exact scattering ampli- 
tude is found to be significantly different with the one in 
Eq. (fT]), then the previous theories should be modified. 
Especially, the separable two-body potential cannot be 
used any more. 

The low-energy p-wave scattering amplitude near 
PMFR has been investigated in Refs. [53, [Hj]. How- 



2 



ever, these studies are based on simplified models of the 
atomic interaction, e.g., zero background potential [53| or 
a separable background potential that decays exponen- 
tially in the momentum space [5o| . The long-range van 
der Waals potential is not taken into account in cither 
case. Due to these simplifications, the scattering ampli- 
tudes given in Refs. [Hj, [EH automatically have the form 
of Eq. JTJ), and cannot be used to judge the applicability 
of the effective-range theory. 

In this paper, based on the realistic long-range inter- 
atomic potential, we provide an explicit calculation for 
the low-energy p-wave scattering amplitude of two spin 
polarized fermonic atoms near a PMFR, and then discuss 
the condition under which the effective-range theory can 
be used as a good approximation. We get sufficient con- 
ditions for the effective-range theory, and show that for 
the ultracold gases of of 6 Li and 40 K with the Fermi tem- 
perature of the order 1/iK, the effective-range theory can 
be used as a good approximation in the resonance regime 
where the scattering volume is enhanced 



A. Main results 

The main results of this paper are summarized as fol- 
lows. 

In this work we first calculate the exact expression 
of the low energy p-wave scattering amplitude with a 
PMFR. We prove that the scattering amplitude can be 
expressed as 



flm,( k ) 



^ + V cff (fe;B;m i )fc 2 + 5 rft \k;B;m z )k + 7^'* (fc;B;m*) 

(2) 

with B the strength of the magnetic field applied along 
the z axis. Here m z is the ^-component of the angu- 
lar momentum; V cS , S cS and lZ cS are fc-dependent scat- 
tering parameters. We obtain the general expressions 
for them. It is pointed out that, the denominator of 
fim, {k) cannot be expressed as a Laurent series with k- 
independent coefficients because V cff , S cS and lZ eS are 
not analytical functions of k. 

Equation @ shows that, in the presence of a PMFR 
the inconsistency of the scattering amplitude with the 
effective-range theory is displayed in a more complicated 
manner. The low-energy p-wave scattering amplitude in 
Eq. @ is different from the one obtained from the stan- 
dard effective-range theory in the following two senses: 

f. The scattering parameters (V eff , iS eff , lZ eS ) depend 
on the incident momentum k. 

2. The term l/[S eS (k; B; m z )k] cannot be included in 
the effective-range theory. 

After obtaining the p-wave scattering amplitude under 
a PMFR, we discuss the applicability of the effective rang 
theory as an approximation of the scattering amplitude 
©. We find that, in the BEC side of the PMFR where 



V eff and TZ eS have the same sign, sufficient conditions for 
the validity of effective-range theory are r^ra « 1. In 
the BCS side of the resonance, the sufficient conditions 
become r\ 1 r2,r^ « 1. Here r\, and are defined in 
Eqs. (|38|) , ([33]) and (|42[) . If these conditions are satisfied, 
the scattering amplitude of our system can be approxi- 
mated as 
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^ + V rff (0:S;m z )fc 2 + TC rft (Ch-B;m,) 



(3) 



which has the same form as Eq. JIJ derived by the 
effective-range theory. 

Qualitatively speaking, the above sufficient conditions 
means that we can use the effective-range theory if the 
fermonic momentum of the cold gas is low enough, the 
magnetic field is tuned close enough to the resonance 
point and the background scattering potential in the open 
channel is far away from the zero-energy shape resonance 
point. For the realistic cold gases of Fermi atoms, if the 
background scattering is far away from the shape reso- 
nance, the effective range theory can usually be used in 
the total region where the p-wave interaction is negligi- 
ble. 

The paper is organized as follows. In Sec. II, we calcu- 
late the p-wave scattering amplitude near a PMFR, and 
obtain the low energy expansion in The parameters 
(V off , S cS , lZ cS ) are expressed in terms of the background 
scattering parameters and the magnetic field. In Sec. 
Ill we discuss the range of applicability of the effective- 
range theory. We show that, under some simple sufficient 
conditions ri,r2 << 1 or ri,r2,f3 << 1, the scattering 
amplitude given by the effective-range theory is a good 
approximation for the exact one obtained in Sec. II. We 
further show that these conditions are well satisfied in 
the cold gases of 40 K and 6 Li when the scattering vol- 
ume is enhanced by a PMFR, and then the previous re- 
sults based on the effective-range theory are applicable 
for these systems. In Sec. IV there are some conclusion 
and discussions. We describe some detail of our calcula- 
tions in the appendixes. 



II. LOW-ENERGY SCATTERING AMPLITUDE 
NEAR THE P-WAVE FESHBACH RESONANCE 

A. p-wave phase shifts with PMFR 

In this section we calculate the p-wave scattering am- 
plitude in presence of a PMFR induced by a magnetic 
field along the z direction. We begin with the two- 
channel Hamiltonian for the relative motion of two atoms 
(Fig. I): 



H 



f+V (b ^(r) 
W(r) 



T 



W(r) 



e(B) 



(4) 
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FIG. 1: (color online) The two channel model of the p-wave 
Feshbach resonance. 



where T is the kinetic energy of relative motion, T/( b &) (r) 
is the background scattering potential in the open chan- 
nel, W(r) is the coupling between the open and closed 
channel, and V^ ol '(f) is the interaction potential in the 
close channel, that has a independent positive threshold 
s(B). In this paper, for simplicity, we assume that the 
background potential V hs (r) is independent on the di- 
rection of r and invariant under the SO(3) rotation. We 
further assume that, in the closed channel there are only 

three bare p-wave bound states l^res*" 1 ) which are near 
resonance with the threshold of the open channel. Here 
m z = 0, ±1 is the angular momentum along the z axis. 

The self energy Em} (B) = /i ros (B - B^ } ) of |^™^) is 
determined by the strength of the magnetic field. The 
difference between BiS*^ with m z = 0,±1 depends on 
the atomic magnetic dipole. For the atoms with small 
magnetic dipole, e.g., 6 Li, the values of Bres and Bits 
are close with each other (53|, while for the the atoms 
with large magnetic dipole, e.g., 40 K, the difference be- 
tween Brel and Brts 1 * 1 is quite large [1]. 

The p-wave scattering amplitude in the open channel 
can be defined with the standard scattering theory [66| . 
To this end we firstly introduce the two-component sta- 
tionary scattering state 




\k) 




where | k) is the eigenstate of the atomic relative momen- 
tum with eigen value hk, and J7+ is the M0ller operator 
defined as 

0+ = Km e-iHT/hjfr/h^ ( 6 ) 

r — >oo 

In the large interatomic distance limit r — > oo, the 
asymptotic behavior of the state l^^) can be expressed 
as 



/ r1 $(+)\ = 1 ( e^+/(f,fc)^ 
V ' k ' (27r»)i V 



(7) 



with r = f/r and f(f, k) the scattering amplitude which 
can be further expanded in terms of different partial 



Ak)Yr*{r)Yr*{ky 



(8) 



Here Y ; m *(f) is the spherical harmonic function. For a 
scattering potential with SO (3) symmetry, the partial 
wave scattering amplitude only depends on the quantum 
number of the angular momentum I. In our case, the 
SO(3) symmetry is broken by the interaction between the 
atomic magnetic dipole and the magnetic field. Then we 
have a m z -dependent scattering amplitude fim z {k). 

In the case of low-energy scattering between two spin 
polarized fermonic atoms, one can neglect all the high 
partial wave scattering amplitudes fim z (k) with I 2, 
and only consider the p-wave amplitudes fim z (k), which 
can be further expressed in terms of the p-wave phase 
shifts 6i rnz (fc): 



Am, 0) 



1 



ik — k cot Sim z (k) 



(9) 



During the scattering process, the bare bound state 
|(/>r«H) is coupled with the p-wave background scattering 
states in the open channel and significantly change the 
p-wave scattering amplitude fim z {k). This effect can be 
directly treated via the Feshbach resonance theory, e.g., 
the methods in Ref. and Ref. 

After a straightforward calculation in Appendix A, we 
find that the final phase shift <5i m> (fc) is the sum of the 
background phase shift S^ s \k) for the background po- 
tential ]/( b s) (r) and a correction A lmz (fc) given by the 
closed channel: 



5 1 (k) = 6<? s) (k) + A lmx (k)- 



(10) 



Here Ai ma (fe) satisfies 



(bg+)\ 12 



(11) 



(5) with ) and 9m z 2 ) given by 



\ mhk 



i4 bg+) ) 



El^lS)^^)*: (12) 



g mz (k 2 ) = Re^lWG^ik^Wl^)- (13) 

In the above we have used the background Green's func- 
tionGV bs) (fc 2 ) 



G^ s) (fc 2 ) 



h 2 k 2 /m + iO+ -T-V (h ^ 



and the background scattering state |0£ bs+ ^} defined as 
|^ bg +>> = \k) + G { + 6 \k 2 )V^\k). (15) 



(14) 
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In the following subsections, we evaluate the low- 
energy expression of the scattering amplitude j\ mz (k) by 
expanding the term — fc cot 8\ mz (fc) in Eq. {9} in the limit 
k — > 0. As shown in Eq. (|10l) . the phase shift Si mz (k) is 
the sum of <J< bg) (jfe) and A l7Uz (k). The low energy behav- 
ior of background phase shift 5^ s \k) is already known 
to be P, |] 



fccot<5< bg) (fc) 



1 1 



1 1 



1 



y(bg) k 2 5 (bg) k ft(bg) 



(16) 



Therefore, if we can further obtain the low-energy ex- 
pansion of term — fccotA lmz (fc), then the expressions of 
— k cot Si mz (fc) and fim z (k) can be calculated straight- 
forwardly. 



is an analytical function of fc, and can be expanded as a 
Talyor series of k: 



F 



(bg) 



Hi 



(r) 



E 



F, 



(bg) 



(2n)! 4-| klm 



(r) 



(22) 



fe=0 



It is pointed out that, all the odd order terms of the 
above Taylor series are exactly zero (66|. Thus the non- 
analytical part of IV^im"^) is included in the term with 
the Jost function ^(k). 

Substituting Eqs. {22]), (03) and {32]) into Eq. {II]), we 
find that in the low-energy limit the factor — fccotAi m _ (k) 
takes the form 



B. The low-energy expansion of — fccot Ai mz (fc) 

In this subsection we investigate the expression of 
— fccotAi mz (k) in the limit k — > 0. To this end, we need 
to expand both the numerator and the denominator of 
(1111) in the low-energy limit. 

In this paper we assume the background scattering vol- 
ume in the open channel is finite. It can be proved that 
(Appendix C), in this case the function g mz (k 2 ) can be 
expanded as 



3^(fc 2 )=5S+^fc 2 + 0(fc 3 ) 



(17) 



(2) 

with g„ lz < 0. 

On the other hand, due to the long-range nature of 
the van der Waals potential, the partial wave scattering 
state IV'ii^?) i s n °t an analytical function of the incident 
momentum k in the neighborhood of k — 66] . To inves- 
tigate the low-energy behavior of IV^im"^) an d then the 
denominator of (jllj) . we separate the non-analytical part 
of IV'fei^?) by introducing the background Jost function 
J{k) ji3| defined as 

Here F^f (r) is the canonical solution of the radial equa- 
tion 



with boundary condition 



p(bg) 
fclm. 



(r -> 0) -> ii(fcr) 



where 



Ji(^) 



cos a; 



(20) 



(21) 



is the first-order regular Riccati-Bessel function |66j . Ac- 
cording to the standard scattering theory 66], F^^ (r) 



fccot A lrriz (k) 



V( A ) (B; fc; m z ) fc 2 ft( A ) (B; fc; m z ) ' 

(23) 



Here we have the fc-dependent parameters: 
V^(B;k;m z ) = -\S{k)\~* ™ m ' - \ 

Mies D-tf 

K^\B;k-m z ) = \jf {k)\~ 2 vw mz x 

" 'h 2 



(24) 



.9 



(2) 



Mros (B — Bq) 



(25) 



with the parameters Bq and w mz defined by 



Bo 



R -o (0) /u • 

"res ilm z l P-res, 

6dfc3 
l d 5 r 

120<ifc5 



(26) 



fc=0 
2 



(27) 



C. The low-energy p-wave scattering amplitude 



In the above subsection we get the expansion (|23|) of 
the factor -fccotAi m2 (fc). Substituting Eq. {23]), (|T6]) 
and (fT0|) into {9J, we finally get the low-energy behavior 
of the p-wave scattering amplitude fi m ,(k) in presence 
of a PMFR: 



ifc 



I 1 I 1 L 1 

~ V a {k;B-m z )k' 2 ~ S ctt (k;B;m z )k ~ Tl ctt (k;B;m z ) 

(28) 
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where the fc-dependent scattering parameters are given 
by 



V eff = 



->eff 



y( b s) 

^(bg) 

y( b e) 2 
1 



1 - 



B-Bq 



v 



iff 2. 



A) 



(1 - 2x + x 2 ) + 



ft( b s) ' 



y(bg 

5^ 



(29) 

(30) 

x 3 ). 
(31) 



not available for 6 Li and 40 K. Nevertheless, the bind- 
ing energies of the p-wave Feshbach molecules are mea- 
sured to be linear functions of the magentic field @, [H| 
in the region with large enough p-wave scattering vol- 
umes (V cS > 10 7 ag). This observation shows that in 
these regions the term \w' mz fires (B — Bq) jw mz | is neg- 
ligible and the effective range R eS can be approximated 
as a fc>-independent constant R^ A \0;k;m z ). 



Here the parameters b m _, and x are defined as 



irw r , 



So far we have obtained the low-energy expression of 
the p-wave scattering amplitude fim z (k) in the case of 
PMFR. With the help of the scattering theory, we ob- 
tain the general expressions (|28M29p for the scattering 
amplitude fim z (k) as well as the scattering parameters 
( V cS , S cS , lZ cS ) , which are formulated in terms of the 
background scattering parameters (V^ hs \ <S' bg ), TZ^^) 
and the magnetic field. It is pointed out that, although 
due to the long-range nature of the van der Waals poten- 
tial we can not express the denominator of f\ mz (k) as an 
Laurent series with fc-independent coefficients, we suc- 
cessfully include all the fc-dependence of the parameters 
(V cS ,S eS ,TZ eB ) into the Jost function J{k). 

Eq. ([29]) clearly show the effect of PMFA that the 
scattering volume V cS diverges under the magnetic field 
B = Bq. In the realistic cold atom gases, for the creation 
of an observable effects with p-wave interaction, the scat- 
tering volume \ V cS \ should be large enough. Particularly, 
|ycff|i/3 should De mu ch larger than the van der Waals 
length so that in the BCS region, the transition tempera- 
ture T c ~ (E F /k B )exp{-TT/(2k%V cS )} [nl of superfluid 
is realizable and in the BEC region the binding energy of 
the p-wave Feshbach molecule be roubst with respect to 
the detail of the atom-atom interaction potential. 

In the end of this section, we consider the depen- 
dence of the effective range R cS on the magnetic field 
B. According to Eq. (f2U|) . R depends on B through 
the ratio x between V^ hgS> and V cS , and the quantatity 
R^(B;k;m z ). In the cold gases of 6 Li and 40 K, the 
background scattering volumes y( bg ) are of the order 
of (10 5 — 10 6 )aQ. According to our above discussion, 
they are too small for the creation of p-wave superflu- 
ids [7l|. Therefore in these systems the strong enough 
p-wave interactions can only be obtained in the reso- 
nance region with \V cii \ » \V^\ or x « 1, which 
implies R cS « R^ A \B]k]m z ). On the other hand, ac- 
cording to Eq. (1251) . the dependence of R^ A '(B;k;m z ) 
on B is significant when the magnetic field is far away 
enough from from the resonant point Bq so that the 
factor Iw^^rcs (B — Bq) /w mz \ is comparable or larger 

than \h 2 / m — gmi\- The values of w' m and w m , are 



III. THE APPLICABILITY OF THE 
EFFECTIVE-RANGE THEORY 

In the above section, we have obtained the expres- 
sion (|2"5)) of the p-wave scattering amplitude fim z (fc) in 
the region near the point of PMFR. It is apparent that, 
this expression is different from the one ([T]) given by the 
effective-range theory in the following two senses: 

1. In the standard effective-range theory, the scat- 
tering volume V and effective-range 1Z are independent 
on the incident momentum k. Nevertheless, in the ex- 
pression (|2"5)l the scattering parameters (V off , S eS , lZ eS ) 
depend on k through the Jost function ^ (fc). 

2. The term 1/[S (k; B; m z )k] cannot be included in 
the effective-range theory. 

It is apparent that, if under some condition the scat- 
tering amplitude (f28f can be approximated as 



1 



ik 



"r V Dtt {0;B;m z )k 2 ^ H ott (0;B;m z ) 



(32) 



i.e., both the fc-dependence of \^{k)\ 2 and the term 
l/[S cS (k; B;m z )k] can be neglected, the behavior of 
the system would be approximately described by the 
effective-range theory. In this section, we investigate the 
conditions for the approximation (|32|) or (??), or the ap- 
plicability of effective-range theory. We will consider the 
importance of the term l/[S cB (fc; B; m z )k] and the fc- 
dependence of the Jost function respectively. 



A. The fc-dependence of the scattering parameters 

In this subsection we search the sufficient condition for 
the ignorance of the fc-dependence of the scattering pa- 
rameters (V eS ,S eS ,TZ eS ). As we have discussed above, 
the fc-dependence of the scattering parameters comes 
from the mode square of the Jost function ^/(fc). In 
the ultracold gases of the fermionic atoms, the maximum 
value of the relative momentum of two atoms is on the 
order of the Fermi momentum kp. Therefore, the impor- 
tance of the fc-dependence of the parameters parameters 
(V , S , 1Z ) can be described by the factor 



ri 



\jjk F )\- 2 -\jm 



(33) 
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FIG. 2: (color online) The functions Df g (k) (blue dash-dotted 
line), Dff(k) (open circle), D g f(k) (black solid line) and 
Dgg(k) (red dashed line) in Eq. (32). 



Obviously, when r% « 1, we can replace \J'{k)\ 2 with 
^(0)|~ 2 and neglect the fc— dependence of the param- 
eters (V cff , S cS , 7e cff ) . 

To investigate the behavior of the the ratio r\ , we first 
calculate the Jost function ^(k). By means of the quan- 
tum defect theory |74j . we can obtain the expression of 
\ a f{k)\~ 2 (see Appendix D): 



k F (1/p 6 ) 



FIG. 3: (color online) Ratio n defined in Eq. (|33[) as a func- 
tion of the Fermi momentum k F with K° =1 = 0.01 (V (bg) = 
-11.79)31, blue solid line), K? =1 = 0.1 (V (bg) = -1.27/9|, 
green empty square), j£Tj° =1 = 1 (V (bg) = -0.23/3f, black 
empty circle), j£Tj° =1 = -0.01 (V (bg) = 11.51/3f, blue dash- 
dotted line), ifj° =1 = -0.1 (V (bg) = 10.4/3|, green filled 
square), K? =1 = -1 (V (bg) = 0, red dashed line). 



{D ff (k) - K? =1 D gf {k)) 2 + (D fg (k) - Ki =1 D gg {k)Y 



where a is a fc-independent coefficient and Dij(k) = 
(k^ 2 Z lj (k) = f,g) with Zij(k) defined in 
In Fig. 2 we plot the functions Dy(fc) in the low-energy 
case. 

The parameter Kf =1 is denoted as Kf =1 is related to 
the background scattering parameters [2| ■ Expanding the 
p-wave phase shift in Eq.(7) of Ref. [2J, we can express 
(V( bg ),S( bg )) in terms of Kf =1 : 



y(bg) = 
5(bg) 



18^° =i r[3/4] 2P6 ' 
35(1 



324(X, =1 ) 2 r[3/4] 4 
The above expression shows that 



Pi 



(35) 
(36) 



Kf =1 is large or the background scattering potential in 
the open channel is far away from the shape resonance, 
then according to Eq. (|54]) . the Jost function \^(k)\ is 
dominated by the term with K^ =1 D g f{k). In this case 
the variation of \ a ^{k)\~ 2 with respect to k is negligi- 
ble and we have r\ « 1. On the other hand, if K? =1 
is close to zero and the background scattering potential 
(34)is close to the shape resonance, then \ c /{k)\ becomes a 
rapid changing function of k and the ratio r\ would be 
significant. 

The above argument is quantitatively verified by Fig. 
2, where the ratio T\ is plotted as functions of k with re- 
spect different values of K® =1 or y( bg ). It is clearly shown 
that, if the fermonic momentum kp ;$ O.l/?^ 1 , then the 
ratio ri and the ^-dependence of the scattering parame- 
ters can be neglected when the background potential is 
far enough from the shape resonance so that y( bg ) < 
If kF ;$ O.Ol/Sg -1 , this restriction can be further relaxed 
to V (bg) < 10/3|. 



B. 



The importance of the term l/[S cfi (k; B; m z )k] 



when Ki_ 



> 



7r/(18r[3/4] 2 ) - 0.1, we have V< bg ) ~ /3| and the back- 
ground scattering potential is far away from the shape 
resonance; when -KjLi is much smaller than 0.1, the back- 
ground potential is in the shape resonance region which 
gives 

Now we consider the features of the ratio r\ , which is 
determined by the parameter K® =1 . Fig. 2 shows that in 
the low energy case with k(3$ << 1, the function D g f{k) 
is almost a fc-independent constant and much larger than 
the other three D-functions. Therefore, if the parameter 



Now we discuss the importance of the term 
1/[S (k;B;m z )k]. Since the purpose of this paper is 
to obtain the sufficient condition for the effective-range 
theory, or the approximation (|32[) . for simplicity, in this 
subsection we assume the condition r± « 1 obtained in 
the above section is already met, and the fc-dependence 
of the coefficients (V eS ,S eS ,TZ eS ) can be neglected. 
In the BEC side of the PMFR where B < B , since 
< 0, the parameters 7Z oS and V cff have the same 
sign. In that case, if the absolute value of l/(<S eff fc) is 
much smaller than the one of l/(V off fc 2 ), then it would 



(2) 
.9 m- 
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be also much smaller than l/(V eff fc 2 ) + l/7Z cS , and can 
be neglected. It is obvious that, in the limit k — > 0, the 
term l/(S cS k) would be much smaller than l/(V cS k 2 ). 
Therefore the importance of the term l/(S cS k) is actu- 
ally determined by the ratio r 2 between the two terms 
for the upper limit of the relative momentum k F : 



'2 



l/[S eS (0;B;m z )k F ] 



l/[V eS (Q;B;m z )k 2 F ] 



(37) 



When r 2 is much smaller than unity, we can neglect the 
term 1/(5 k) . If r 2 is comparable or larger than unity, 
the term 1/(5 k) would be necessary for the theory. 
The straightforward calculation with Eqs. (|2T))) . (13TH) . 
and ((Ml) yields 



''2 



0Z 



35\V° s (k F ;B;m z )\ 



(p 6 k F ). 



(38) 



In the practical cold atom systems we have f3§k F « 1. 
Therefore, in the resonance region with V off > /9|, the 
ratio r 2 in Eq. Q38p is much smaller than unity, and then 
the term 1/(5 k), is negligible. 

In the BEC side of the resonance with B > Bq, the 
parameters V and 1Z have different signs. In that 
case, there is a speical momentum 



11^(0; B;m z ) 
V cS {0;B;m z ) 



(39) 



which makes the terms 1/[V k 2 } and 1/[1Z ] cancel with 
each other or 



1 



1 



V cS (0;B;m z )k 2 K eS (0;B;m z 



0. (40) 



Therefore, if the atomic relative momentum k is far 
away from the absolute value of l/(1Z cS k 2 ) would be 
quite different with the one of 1/V eff . In that case we 
can still neglect the term l/(5 cff /c) under the condition 
r 2 « 1 or |l/(5 off fc)| « \l/(Jl cfi k 2 )\. 

If the atomic relative momentum k is in the neigh- 
borhood of and the terms \/{lZ cS k 2 ) is canceled with 
1/V eff , the scattering amplitude (1251) can be expressed as 



1 



S c « (0;B;m z )k, 



(41) 



In that case, if the absolute value of 1/ [5 eff (0; B; m z ) £;*] 
is much smaller than k*, we an also neglect the term 
with 5 eff , even in the neighborhood of fc*. We define a 
parameter r% as 



''3 



|5 eff (0;B;m z )k 2 



Then the term with 5 can be neglected when r 2> 3 << 1. 
A further calculation with Eqs. (j39)) . ([2"9 ]) . (j30)) . f?5l) and 
(|36| implies that 



7T/3| 



35|V off (0; B; m z ) TZ eS (0; B; m z ) \ ' 



(42) 



As shown above, the condition r$ << 1 is obtained 
for the momentum region k ~ fc*. Since the realistic 
momentum of the atomic relative motion takes the value 
between zero and k Fl in the cases with k F < k*, we can 
disregard the restriction of the ratio , and use effective- 
range theory under the condition rj.,2 << 1 in both the 
BEC and the BCS sides of the resonance. 



C. Summary 

In summary, the general sufficient conditions for the 
effective-range theory in the BCS side of the resonance 
can be summarized as 



ri,r 2 ,r 3 « 1, 
while the ones for the BEC side are 
ri,r 2 << 1. 



(43) 



(44) 



From the definition of the ratios n , r 2 and , we no- 
tice that in the realistatic cold gases of Fermi atoms, the 
crucial factors for the usage of effective range theory is 
the background p-wave scattering volume V^ bg ^ and the 
B-dependence of the factor (0;B;m z ). If the the 

background p-wave scattering is far away from the shape 
resonance so that V^ 8 -* ~ /3| ~ (100a ) 3 then according 
to our previous discussions and Eqs. (|38|) . the condi- 
tions n , r% « 1 can be satisfied in the region V > 
10/3| ~ 10 7 aQ where the p-wave interaction is strong 
enough for the creation of p-wave superfluids. In that 
region we also have lZ eS (0; B; m z ) — (0; B; m z ). 

If {0;B;m z ) can be further approximated as a B- 
independent constant which is of the order then the 
condition r% << 1 can also be satisfied, and the effective 
range theory can be used a good approximation for the 
realistic scattering amplitude. In the following subsec- 
tion we show that the PMFRs in the cold gases of 40 K 
and 6 Li are just of this case. 



D. Discussion for the cold gases of K and Li 

In the above subsections we obtained the sufficient 
conditions (T4"3"|) , (|4"4"1) of the effective-range theory for 
the p-wave scattering amplitudes of polarized fermonic 
atoms near a PMFR. In this subsection, with the help of 
the conditions, we perform a discussion on the usage of 
effective-range theory in the ultracold gases of 40 K and 
6 Li. 

For the ultracold gas with 40 K atoms in the state 
|9/2, —7/2), we have C 6 = 3897(a.u.) [z| and V( bg > = 
— 10 6 ao 5]. These parameters leads /3q = 130ao and 
Ki =1 = —0.16. If the Fermi temperature T F = 1/xK, 
then we have k F f3s = 0.06. The straightforward calcu- 
lation shows that rj = 0.01. Therefore fc-dependence of 
the scattering parameters can be safely neglected. The 
p-wave Feshbach resonance for the states with mi = ±1 
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occurs at B = 198.37G with width AB = 25G and 
effective-range lZ cS — 47.2ao- The resonance for the 
states with mi — occurs at Bq — 198. 85G with width 
AB = 22G and effective-range K cS = 46.2a [5]. Ac- 
cording to these data we have < 0.02 when fe* < kp. 
Then effect from the ratio r% is also negligible. Therefore 
the sufficient condition for the usage of effective-range 
simply becomes r 2 « 1 . Further calculation shows that 
r 2 < 0.013 when |V cff | > |V( bg )|. Then the effective- 
range approximation (I3"2l is applicable for 40 K atoms 
in the state |9/2, -7/2) in the whole region of PMFR 
with |V off | > |V (bg) |. The condition for the effective- 
range approximation is broken only in the small region 
220.5G < B < 221G (mi = 0) or 223G < B < 223.7G 
(mi = ±1) where we have |V eff | < 0.005/3| or r 2 > 1. 

Now we consider the gas with 6 Li atoms in the ground 
hyperfine state \F = l;mp = 1). In that case we have 
G 6 = 1393(a.u.) [zSlzl and V( bg ) = -(35.3a ) 3 . These 
parameters leads /3q = 62ao and Kf =1 = —0.38. If the 
Fermi temperature Tp = 1/iK, we have kpfte = 0.01 
which implies r± — 2 x 10~ 4 . Then similar as above, the 
effective range approximation (|32l) is also applicable for 
Li atoms in the whole region of PMFR with |V cff | > 
|V (bg) |. 

IV. CONCLUSION 

In this work we obtain the explicit expression of the 
p-wave scattering amplitude of two ultracold spin polar- 
ized fermonic atoms near the p-wave Feshbach resonance. 
We show that due to the long rang nature of the van 



der Waals potential, the scattering amplitude is explic- 
itly described by Eq. (2J in the low-energy case. With 
the help of the quantum defect theory, we formulate all 
the scattering parameters (V eff , <S eff , lZ eS ) in terms of the 
background parameters and the inter-channel coupling. 

Based on this result, we discussed the applicability of 
the effective-range theory, or the Eq. ((3]) as an approxi- 
mation of the exact scattering amplitude. We show that, 
in the BEC side of the resonance, the sufficient condi- 
tions of the effective-range theory can be quantitatively 
described as ri, r-x « 1 while in the BCS side the condi- 
tions become T\ , r<x , << 1 , where r\ , r% , r 3 are defined 
in Eqs. (35]), (33J), and (22). The applicability of the 
effective-range theory for the ultracold gases of 40 K and 
6 Li arc examined with our results. The effective-range 
theory is shown to be a good approximation in both of 
the two cases in the resonance regime where the absolute 
value of the scattering volume is equal to or larger then 
the background one. 



Appendix A: The p-wave phase shift with PMFR 

In this appendix we post the derivation of the p-wave 
phase shift 8\ mz (k) in Eqs. (fTUf and (1TT1). Our calcu- 
lation is based on the method in Ref. [60|. We begin 
from the scattering state ) in Eq. (5]). According to 
the scattering theory (66j . the open channel component 
^(op+)^ an( j c i osec i channel component |0g Cl ' 1 ) of |$g .) 
satisfy the two channel Lippmman-Schwinger equation 



i4 op+) > 
i4 cl+) ) 



\k) + G^ +) (fc 2 )0 bg >|4 op+) ) + G( +) (fc 2 )IF|4 cl+) ); 

4 cl) (fc 2 )^i4 cl+) )+4 cl) (fc 2 )PFi4 op+) ) 



(Al) 
(A2) 



with the free Green's functions 



4 +) (fc 2 ) = 

G^(fc 2 ) = 



1 



h 2 k 2 /m + iQ + — T 
1 

h 2 k 2 /m-f-e(B)' 



(A3) 
(A4) 



We further define the background Green's function 
G^)(k 2 ) and the closed channel Green's function G c \ as 



(A5) 



H 2 k 2 /m + i0+ -f -y( bg )' 

G c i(fc 2 ) = ^ . (A6) 

h 2 k 2 /m-T -y( cl ) -e(B) 



Then we have the relationships 

G^\k 2 ) = G[ ) + ) (fc 2 )-G[ ) + ) (fc 2 )^ bg )GW(fc 2 )(A7) 
G^(k 2 ) = G cl (k 2 )~G d (k 2 )V^Gt\k 2 ). (A8) 



Substituting Eqs. (IA7|) and (|A8|) into the last terms of 
the right hand side (r.h.s.) of Eqs. (|A1[) and (|A2|) . and 
using the Lippmman Schwinger equation (j!5[) for the 

background scattering state \<jf~ s+ ^), we get the equa- 
tion which relates ) with the background scattering 
state |0? g+) ) HI: 



i4 op+) ) 



i4 bsH 



) ) + Gi + >(k 2 )W\^) (A9) 



14 ) 



G cl (fc 2 )^|4 op+) ) 



(A10) 
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To calculate the p-wave phase shifts Si m ^(k), we 
operate the projection operator V m , for the manifold 
(I = 1, L z = m z ) on both of the two sides of Eqs. (|A9j) 
and (|A10p . Then we have 

Kul4 op+) ) = ^4 bg+) > + Ku<V)w^ m j4 cl+) ) 

(All) 

r m M ) = <t>™ z) >— < \ ■ ( A12 ) 

Here we have used 

PtiuGfiV) = Pm.GffC* 2 )^,; (A13) 
Pm z Gci(fc 2 ) = P mi Gd(fc 2 )P mi (A14) 

which are guaranteed by the rotational symmetry along 
the z-axis of the system. We also made the approxima- 
tion 



H 2 k 2 /m - fi ICS (B ~ Bt z) ) 



(A15) 



That is, we only take into account the contribution from 
the near- resonance bound state \4>r^ z ^) in the closed 
channel. 



Substituting Eqs. (|A12[) into Eq. (|A11I) . we get 



•P,,.l4 OP+) > = 7>m.\<t>$* +) ) 



+G i £(k 2 )W\4> ( i Z* ) )A {m > ) (B,k 2 ) 



with 



A^\B,k 2 ) = 



d m *) 



h 2 k 2 /m-n iea (B-B^) 



(A16) 



(A17) 



Replacing the V m J<f>f p >) in the r.h.s. of Eq. (jA"T7|> 
with the r.h.s. of (IA16|1 . we get 

A^{B,k 2 ) 

_ (^ ] \W\^ +) ) 

~ h 2 k 2 /m - ^ ICS (B - Bt z) ) - {4T/ ) \WG { h ^(E)W\4Z' ) ) ' 

(A18) 



Substituting Eq. (|AT8| into Eq. ([176]) . and using the 
asymptotic expression J7| of the scattering state and the 
definition (|5J| of the partial wave scattering amplitude, 
we can obtain the p-wave scattering amplitude 



J 



* h 2 k 2 /m - /i ros (B - Bt* ] ) - (4&l M) \WG£ ) (k 2 )W\<l& 

r 



ft M - f (bS) (k) - - Whim, lyres /xyrca |^ l^lm, / ( > q] 



Here l^ii^) is defined in Eq. (H). IV^imj) is defined as 
as 



^K^)^ 1 ^^ (A20) Wfe) = /lW(fc) -^ W iiM^ (A23) 



with 



4 b s-)) = |fc) + i — , --^ (bs) |fc) 

k h 2 k 2 /m + i0- -T -W b s) 



with 



(A21) C (k) = * x 



the background state with in-going boundary condition. (22 ( ) 2 \ 

In the above calculation we also used the asymptotic be- ( /i fc /m — /i ros 

havior of the background Green's function: 

lim (fiG^ik 2 )^) = -mS^-- (^i b - g_) |r) (A22) In the derivation of Eq. (|A23[) we have used the relation- 

r ^°° S \ 2h r r ship (Appendix B) 

with f — r/r. 

With straightforward calculation, we can further 
rewrite the scattering amplitude fi mz [k) in Eq. (|Al9|l HAkim ) = e~ 2l<5lbg > ( A25 ) 
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and 

G { h V(k 2 ) 



(mh) / dk 



, r J^ g+) )(^ S+) | 
k 2 + i0+~k> 2 

2 



= -(m 2 /h)iri I dk'5{k 2 -fc' 2 )|^ g+; }(<^ I 



(A26) 



Here P J refers to the principle value of the integral. 

Considering the relationship ((9]) between the scattering 
amplitude fim z (k) and the phase shift 5i mz (k), it is easy 
to prove that the phase shift 5i m ^(k) corresponding to 
the scattering amplitude (|A23I) is the one given in Eqs. 
(|TU| and (fTTj) . It is pointed out that, this result can also 
be proved with the method in Ref. f[67|V 



Appendix B: The scattering states with ingoing and 
outgoing boundary conditions 



In this appendix we prove the formula (|A25|) in Ap- 
pendix A. We begin from the relationship [69j between 
the three dimensional scattering states with ingoing and 
outgoing boundary conditions: 



(Bl) 



Considering the definitions (H]) and (TAlOl) of 
we can obtain 

E^K^'C*)* = Y,(A~M Y r°(-k) (B2) 



We further define the one dimensional functions 
i^ bg±) (r) as [66| 



Using the relationships 

Y^'i-k) = (-l) m * +1 yf m *(fc)*; 
Y^'{f) = (-l) m *yf m ' (r)* . 



we get 



(bg-). 



(B3) 



(B4) 
(B5) 



(B6) 



On the other hand, we know that F fc (bg+) (r) and 
F ki 8 ~ ] ( r ) satisfy the same differential equation 



with the same boundary condition F^ s±) (0) = 0. Then 

F ki S '\ r ) is proportional to F^ s+) (r). To calculate the 
ratio between FS^^ (r), we consider their asymptotic 
behaviors in the limit r — > oo: 



(B8) 



^iTV) = h(kr) + kf^'ikYlMkr) - iji(fcr)] 



(bg). 



with ji(a;)defined in ([2"T]) and 



„ , . cos a; 

n\{x) = sini 



(B9) 



(BIO) 



the irregular first order Riccati-Bessel function. Compar- 
ing the coefficients of n;(fcr), we have 

F^-\r) = e-wfl^-lfr). (Bll) 
Here we have used 



f^\k) 



ik -k cot <^ bg) (k) 



isin^ bs) (fc)e< g, ( fe ). (B12) 



Substituting Eq. (|B11|) into Eq. (JESJ , we get Eq. (|A25|) . 



Appendix C: The expansion of the factor gm^{k 2 ) 

In this appendix we prove the Eq. (fl7j). We firstly 
rewrite the factor j mz (fc 2 ) as 



9mz (k 2 ) = Re(^\WG^\Q)W\^ ] ) ~ fc 2 Re(^WGi bg) (0)GV bg) (fc 2 )H/|^)>. 



(CI) 



Here we have used the identity 

G<L bg) (fc 2 ) = G^ bg) (0) - k 2 G^ B) (0)G^ s \k 2 ). (C2) 



The first term in the r.h.s. of (|Clj) is independent on k 
It contributes the constant term gml in Eq. (TTTl) . 
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On the other hand, the second term in Eq. (IC1[) can 
be re-written as 



fc 2 Re( 



'- ) \WG { + e \0)G i + e) (k 2 )W\(l)(™* ) } 

ft 4 j™0+ J (fc' 2 - - k 2 - k 2 ) ■ 



(C3) 



In the limit k = 0, we have 

lim Re(4Z x) \WG^ e) (0)GV bg) (fc 2 )W^ 



dp 



(CM) 



We know that the function K^S^I^I^i^)! 2 decays 
to zero when p — > oo. On the other hand, as we have 
shown in Sec. III.B, the factor \ a ?{k)\ 2 tends to a non- 
zero constant in the low energy limit if the background 
scattering volume in the open channel is finite. Using the 
relationship ggj between J(k), i^ (bg±) (r) and I^S?)' 
and the low-energy behavior dHJ) of^ gi) (r), it is easy 
to prove that the factor |(</>^ z ' ) |MK|'i/'pi^)| 2 is propor- 
tional to p 3 in the limit p — >• 0. Therefore, the above 
integration in Eq. (jC4|) converges to a finite constant in 
the limit k 2 — > 0. Then the expansion in Eq. (JTTJ) is 
proved and we have 

9%1 = -(4^ ) \WG^\0)G^ S) (0)W\^ ) ) < 0. (C5) 



(,(bg+)\|2 



Appendix D: The background Jost function 

In this appendix we calculate the Jost function ^ (k) 
of the background scattering state. To this end, we 



introduce a function F, 



(bg±) 



(r) 

F ki S±) ( r ) is defined in Eq. (T 



F^ s± \r)/k 2 where 



It is apparent that 



-^fci S± ^ ( r ) ^ s a solution of the radial equation 



dr 2 



y(bg)( r ) 



F^\r)=k 2 F^\r) (Dl) 



(bg). 



with a fc-independent boundary condition 



(D2) 



Following the spirit of quantum defect theory [74[ , we 
assume the scattering potential V^ bg ) (r) can be approxi- 
mated as the van der Waals potential — ft 2 /3g/(r 6 m)when 
r is larger than a critical distance ro which is much 
smaller than /3g. In the region with r < ro, V^ hs '(r) is as- 
sumed to be so large that the atomic kinetic energy fc 2 is 



negligible in comparing with V^ bg )(r), and then F^ s \ 



is independent on k. 



In the region r > ro, F^ (r) is the superposition of 
the two independent solutions Xei ( r Q an d K fi ( r ) °f Eq. 

(in Ref. [68j], Xei' 1 ( r ) an d i^ e \ (r) are denoted as 



jftr) and ^ } (r)): 



^ SJ (r) = a fcX ^ 1 (r)+/3 fc ^ 1 (r) 



(0) 



(0) 



(D3) 



In the short distance region with r << /3q, Xi,3i( r ) an d 
w i.2i( r ) are almost independent on fc [74j |. 

The wave function F^ s ^ (r) in the two regions are con- 
nected at point r = ro, where F^ 6 \r), xii v") and 
(r) are approximately independent on fc. Then we 
know that and /3fc are independent on fc and we have 



with efe = h 2 k 2 /m. The above equation yields 

Fi bg) M = fc 2 axS(r) + fc 2 /3^(r). (D5) 

In the region r — > oo, the asymptotic behaviors of 
Xe?i( r ) and t^i(?*) can be expressed in terms of Zij(k) 
(i, 3 = f,9) defined in [HI 



,(°) 



.(0) 



x£i(r) 



9 cos ( fcr - - 



(D6) 



Z s/ sin (kr - |) - ^gg cos (fcr - |j 



(D7) 



On the other hand we know that in the same limit we 
have 166] 



F^\r) = f{k) [ e < E) ( fe )sin(fcr-| + ^ bs) (fc)) 



(D8) 



Together with Eqs. ([D5|) . (|D6]) . (|D7|) as well as (|D8|) . we 
obtain the expression of (f3"4"|) of |^(fc)| 2 and the back- 
ground phase shift: 



tan 5{ 



(bg) 



Kl=iZga - z fg 



KUZ g} Z 



ff 



(D9) 



where D tj (k) = (fc/3 6 ) 3 / 2 Z„ (fc) and Kf =l = -j3/a is the 
one in Sec. III. The result (|D9[) is also given in Q. 
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